CHIRAL MESON MASSES AT FINITE TEMPERATURE AND 

DENSITY 



A. Barducci, R. Casalbuoni, G. Pettini 

Dipartimento di Fisica, Universitd di Firenze, 1-50125 Firenze, Italy 
and Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, 1-50125 Firenze, Italy 

o" 
o 
o 

CS| ■ R- Gatto 

Oh' ^ 

D ■ Departement de Physique Theorique, Universite de Geneve, CII-1211 Geneve 4, Switzerland 

(N . UGVA-DPT 2000/08-1087 

^ . Firenze Preprint - DFF - 359/06/2000 

in 

oo 

0\ ; Abstract 

O 
O 
O 

Q-i! The ratio of the sigma mass to the pion mass at finite temperatures and 

D , densities provides for a quantitative signal of chiral symmetry breaking. We 



X 



calculate this ratio by using an extension to finite chemical potential of the 
field theoretic composite operator formalism as applied to QCD. The calcu- 
lation is limited to regions of the phase diagram where only quark-antiquark 
condensates dominate (no quark-quark condensates) and it confirms the ex- 
pected behaviours. In particular the sigma becomes an essentially stable 
particle in a narrow region bordering the transition line from broken to re- 
stored chirality . This pattern is qualitatively the same both for the region 
where the transition is of second order as for the region where it is of first 

order (apart from the discontinuities expected in the latter case). 
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I. INTRODUCTION 

Since the original discussion by Nambu [^ in 1960, chiral symmetry breaking has played 
a crucial role in the theory of strong interactions. In 1974 Lee and Wick proposed that 
at high densities chiral symmetry may be restored. Restoration at high temperatures was 
discussed in the same period by Dolan and Jackiw |^ and by Weinberg |^. Condensation of 
quark-antiquark pairs plays the main role in chiral symmetry breaking, except for very high 
densities where quark-quark condensates are expected to play the important role. In this 
note we do not consider this very high density region and limit ourselves to densities such 
that only quark-antiquark condensates are important. Such a region is the one important 
in hadronic physics, and chiral restoration from this region is expected to be reachable with 
high energy heavy ion collisions ||^. Chiral transition through the quark-antiquark order 
parameter is also expected to be of great relevance for the physics of the early universe [^. 

A great number of studies has been carried out on the expected QCD phase diagram 
(see for instance 0). For a review of the subject see the recent paper by Rajagopal P] 

Lattice calculations are rapidly progressing and they will with time provide for a complete 
approach (a review can be found for instance in Ref. 0). Many theoretical approaches use 
the Nambu- Jona-Lasinio model and different dynamical approximations for quantitative 
indications on the chiral transition (see for instance ||1CI||). 

Besides the chiral transition also the deconfinement transition takes place. A single 
transition is generally expected to occur, rather than separate transitions for deconfinement 
and chiral-symmetry (see for instance for a heuristic argument Ref lll|). Theoretically one 
has to deal with different order parameters. Moreover they cover extreme ranges in opposite 
domains. The transition from the confined to the non-confined phase is often described in 
terms of the thermally averaged Polyakov loop which applies particularly in the limit of 
infinite quark masses. For the chiral symmetry transition the typical order parameters are 
the thermally averaged quark-antiquark bilinears (at least up to not very high densities) 
which apply mostly in the limit of vanishing quark masses. 



The current quark mass violates explicitly the chiral symmetry whose restoration char- 
acterizes the phase transition. Its role is formally analogous to that of the external magnetic 
field in the ferromagnetic transition. As long as the mass is small, such as for u and d 
quarks, the notion of phase transition can be approximately useful. 

Our study uses a composite operator formalism for both finite temperature and finite 



density. It is based on constructing the effective action for the composite operators |12|jl3[ 
in the presence of temperature and chemical potential. We focus on the study of the ratio of 
the scalar to pseudoscalar mass Mcr/Mj^ in the plane of temperature and chemical potential. 
The ratio might have a physical interest as a possible experimental signature for the tran- 
sition, although it is too early for assessing its usefulness in these regards. It has anyway a 
theoretical interest as it constitutes a significant parameter for the degree of chiral symme- 
try breaking at finite temperatures and densities. The ratio is equal to unity when chiral 
symmetry is restored (we neglect for the moment the small u and d current quark masses). 
In the hadronic phase, at zero temperature and small density, in the broken chirality phase, 
the pion is a goldstone of vanishing mass, while the sigma has a finite mass from the chiral 
condensation. The ratio of the scalar to the pseudoscalar mass is thus expected to decrease 
from very large values to a value of one at the chiral transition (the u and d current quark 
masses only slightly change this picture). 

A physical consequence of this behaviour is that, before approaching the critical values 
for the transition, by increasing together or separately temperature and density, for the 
sigma particle there is not any phase space left for its decay into two pions. Its main source 
of instability is thus suppressed when approaching the phase border in the two dimensional 
plane of temperature and chemical potential. Consequently the sigma becomes an almost 
stable particle with dominant decay channel |T^ that into two photons. The suggested 
experimental signature would thus be in that region a bump in the invariant mass of two 
photons. We do not know at this stage how this will be visible in heavy ion collisions. 

On the other hand it is worthwhile to recall that indications on the sigma at zero tem- 
perature and density have been rather elusive just because of the very large width the sigma 



has in that case. It has not been demonstrated so far that the pion-pion phase shift passes 
through 90 degrees in a region below some 700-800 MeV, suggesting that the large width 
may mask the overall effect. Also it is well-known that strong rescattering effects in differ- 
ent channels make the interpretation of the low energy phenomenology in the O"*""*" channel 
not easy |TB|. We just mention that a recent experiment by the E91 collaboration brings 



evidence for a light and broad scalar resonance in the D-decay channel into three pions ||T^ 



and that a simple model makes plausible the interpretation in terms of a sigma resonance 



17| . For further phenomenological discussions on the sigma resonance in hadronic physics 
we refer to the papers posted in relation to the recent conference on "Possible existence of 
a light sigma resonance and its implications to hadron physics", held in Kyoto, June 11-14 
(for a review see the summary talk by Tornqvist [jl8[). We also mention that there is been 
recently discussion as to the possibility of a glueball interpretation of the sigma resonance 



and for this we refer to a recent talk by Pennington disfavoring such an interpretation [jT9 

We recall that the QCD phase diagram, even for massless quarks, is known to have a 
richer structure when studied in the plane of T and chemical potential, as one expects both 
types of transitions, of second order for low chemical potentials and of first order beyond 
some critical value of the chemical potential. Besides, for large chemical potential a color 
superconductivity phase should occur, but here we limit ourselves to consequences of quark- 
antiquark condensations and assume not yet to be in a region where quark-quark condensates 
appear. The separation along the critical line between the two types of transitions defines 
a tricritical point which has been the subject of extensive studies [^,^. As far as the 



particular subproduct of the present calculation, that is finding the region in the phase 
diagram where the sigma becomes of such a low mass to be no longer able to decay into 
two pions and thus appearing as a narrow state, we find that, independently on the type 
of transition, this region is anyway confined to a small strip near the transition line. What 
changes of course is that when passing through the transition line one finds for the the ratio 
of the sigma to the pion mass a discontinuity in the case when the transition is of first order. 
In Section II we recall the main steps to evaluate the effective potential in the composite 



operator formalism. Then we derive the expressions for the pion and cr-meson masses at 
finite temperature and density in our model. Finally we obtain a formula for their ratio and 
we show the numerical results for typical behaviours when crossing the lines of second and 
first order phase transitions. 

II. EFFECTIVE ACTION 



The zero temperature Euclidean effective action for an SU{N) QCD-like gauge theory 
is (we follow [ |T^ ) 
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where 



So ^ = {ip- m) 
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m is the bare quark mass matrix and r2(S) is the sum of all the two-particle irreducible 
vacuum diagrams with fermionic propagator S and S = —6r2/6S. One has at two-loop 
level 



12 = -Tr(SAS) 
where A is the gauge boson propagator, so that 
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One can therefore rewrite Eq. ([l|) 
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A parametrization for 5] employed in [Q was 



S = (s -f i75P)/(^) = ^s + nsSy 



(6) 



with s and p scalar and pseudoscalar constant fields respectively and with a suitable Ansatz 
for fik). 



From Eq. (||) (see Ref. [|T^]) one obtains the effective potential 
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where C2 is the quadratic Casimir of the fermion representation (for SU{3)c C2 = 4/3), g is 
the gauge couphng constant and Q the four-dimensional volume. Furthermore 
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with a = 0, ■ ■ ■ , 8; Aq = a/2/3 / and Aj (i = 1, ■ ■ ■ , 8) are the Gell-Mann matrices. 
For f{k) we have chosen as a variational Ansatz [^] 



(9) 



fc2 + M2 

where M is a momentum scale expected to be of the order of Aqcd- 

The extension to non zero temperature T and non zero chemical potential /i corrects the 
effective potential for a quark of mass m to the form 
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where k'^ = {ko + i/i, k) and we have used the Poisson formula, equivalent to the application 
of standard imaginary-times rules [Q, p0| , |22 |. 



In Eq. (Iig 



sfCk) 



(11) 



and analogously for Sp. The coupling constant depends on T and /i as 



-^^c(T,^)^co + c,(T,/.) = co + yln(l+e^ + C^^ 



(12) 



with b = 2A7t'^ / {UN — 2nf) , and we will discuss the determination of the parameters cq, M, ^ 
and ( later on. Also for non vanishing T and fi we impose the normalization 
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The following relation relates the minimum of the potential s{T,fi) to the condensates 

N M^ 
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When using physical pion fields with their proper normalization, the normalization con 
dition of Eq. (|13D is equivalent to the Adler-Dashen formula 
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III. MASSES 



To derive the masses we use canonical fields 



and we obtain 
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where s is the minimum of the effective potential in presence of the quark bare mass m. 



In the following we use the minimum condition 
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where S^ = sf{k). Furthermore we have the normalization condition for the effective 
potential 
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with Sq = So/(/c) and Sq is the minimum for the case of a massless quark. The two 
preceding equations allow for the following expressions for the masses 
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The parameters cq and so(T = /i = 0) are fixed by following the procedure of Ref. 

I3|p0| , |23[] where one was able to fix the value Cq = 0.554 from the renormalization con- 
dition and the gap equation for the effective potential in the limit of small current quark 
masses. By inserting this value for cq in the gap equation for massless quarks, one finds 
that chiral symmetry is spontaneously broken (at T = /i = 0). Actually the origin of V^ is a 
local maximum whereas the absolute minimum is at so(T = /i = 0) = —4.06. To determine 
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the mass scale M and the quark masses from the experimental data, one has to derive the 
explicit expressions for the masses and for the decay constants of the pseudoscalar octet 
mesons which are the pseudo-Goldstone bosons of chiral symmetry breaking. These expres- 
sions constitute a system of coupled equations which we have solved by an approximation 
method. The experimental inputs are the decay constant and mass of the charged pion, 
the charged kaon mass and the electromagnetic mass difference between the neutral and 
the charged kaon. The outputs of the numerical fit for the octet meson masses (agreement 
within 3%) are the masses of the u^d and s quarks and the mass parameter M. The values 



we get are the following 1 12,13,23 



M = 280 MeV 

ruu = 8 MeV 

ma = 11 MeV (23) 

rus = 181 MeV 

m = = 9.5 MeV 

2 

Finally, the parameter ^ in the expression ([T2|) for the running coupling constant has been 
set to ^ ^ 1 [24] by comparing our model in the low-T and /x = regime with the results of 



Ref. [25 1, whereas ( is left as undetermined [20,26[. 



With these values for the relevant parameters it is easy to determine the masses of the 
pseudoscalar and scalar mesons at T = /i = (we are using the value /^ = 93 MeV for the 
pion decay constant) 

M^ = 135 MeV 

M„ = 630 MeV (24) 

and (again at T = /x = 0) the fermionic condensates at M for massless and massive quarks 

{i/jiIj) ^ (-197 MeVf 

(iPij) ^ (-200 MeVf (25) 

9 



respectively. 

IV. THE RATIO M^/Ml 

From Eq. ( ^21) we can calculate the ratio M^/M^. We then use the expression for f^M^ 
obtained from Eq. (0), and get 
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As we have said, it is particularly interesting to evaluate when, on the phase diagram, the 
mass of the sigma becomes smaller than twice the pion mass, so that it cannot kinematically 
decay into its main decay channel. The phase diagram is now of a more complicated structure 
than in absence of chemical potential. As we know, we expect a tricritical point (/i(,Tf) 
separating the region of second order transition from that of first order. We thus have to 
consider different cases. For /i > /ij, by letting T grow , one passes through the first order 
branch of the critical line. For yU < /^t, increasing T one passes through the second order 
branch. 

We recall that of the two parameters present in the coupling, C, is left undetermined and 
can be only constrained within some limited region from qualitative considerations. These 
are constraints on the critical values of the chemical potential at very low temperatures, 
which suggest to locate this parameter in the range ^ ^ 0.1 — 0.3 to get reasonable values 
for the critical quark densities ||2^ . At the same time the location of the tricritical point as 
well as the appearance of a minimum in the origin (namely the line of second order phase 
transitions and its continuation at low temperatures) are very few affected by the choice 
of this parameter. To be specific, the coordinates of the tricritical point {nt,Tt) (in MeV) 
change from (76, 92) for C = 0.1 to (75, 91) for C = 0.3. 

10 



To evaluate M'^/M'^ in Eq. (^) we can perform the integration in k^ by using the Cauchy 
theorem and summing up the series; then, through a numerical calculation we can finally 
perform the integration left. Typical behaviours of the ratio M^jM.,^ at fixed /i and growing 
T for fi < fit and fi > fit are shown in Fig. 1 and Fig. 2. We have chosen ( = 0.1, but 
almost the same pictures show for other values of ( in the range 0.1 — 0.3. Thus Fig. 1 is for 
/i = 70MeV{< fit), whereas Fig. 2 is for fi = 120MeV"(> fit). The current quark mass is set 
to m = {rriu +md)/2 = 9.5MeV. The little jump visible in Fig. 2, which occurs after the a 
mass has decreased below twice the pion mass, is due to the jump in the fermion condensate 
when crossing the line of first order phase transitions. The size of this jump increases for 
higher values of the chemical potential, namely by going far from the tricritical point. 

V. CONCLUSIONS 

Clear-cut experimental signatures for the sigma meson in low energy QCD phenomenol- 
ogy are made difficult by its expected large width and by complicated unitarity contributions 
and related distortions. We have examined here the picture emerging at finite temperatures 
and densities. In the simplest model with massless u and d quarks the pion and the sigma 
are degenerate in mass only as long as chiral symmetry holds. At zero temperature and 
in absence of chemical potential, chiral symmetry is broken, the pion is a goldstone of 
vanishing mass, and the sigma has has a finite mass from the chiral condensate. In the 
temperature-density phase diagram one expects that, when approaching, from the chirally 
broken phase, the line of chiral restoration, the pion mass increases and the sigma mass 
decreases, and finally the two particles become mass degenerate after passing into the region 
of chiral restoration. Before reaching the transition line, within a certain strip touching 
such a line, the sigma will not have phase space left to decay into two pions. The dominant 
decay channel will be that into two photons, and in such ideal conditions the signal from a 
sigma will be a narrow bump in the two-photon invariant mass. The ideal picture will be 
quantitatively changed because of the finite u and d quark masses and the role of the strange 
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quark, but the qualitative feature of a narrow sigma resonance is expected to remain there 
even after inclusion of these effects. It is too early to say whether a possible experimental 
signature for the chiral restoring transition may be inferred from such behaviour. On the 
other hand determination of the narrow two-photon bump would be of great theoretical 
interest within the chiral picture. 

The present calculation relies on the composite operator formalism applied to QCD at 
finite temperature and chemical potential. It uses the same QCD model parameters that 
had been derived in previous works, using the same formalism, for the masses and decay 
constants of the pseudoscalar mesons. The experimental inputs were the mass and decay 
constant of the pion, the charged kaon mass and the electromagnetic mass difference between 
the neutral and the charged kaon. The effective potential is consistently renormalized to 
derive expressions for the masses of the scalar and pseudoscalar mesons at finite temperature 
and chemical potential. We have limited the calculation to the phase region where quark- 
antiquark condensates dominate (no quark-quark condensates). The ratio of the sigma mass 
to the pion mass is then calculated within the composite operator formalism for finite tem- 
peratures and densities. We recall that the same model employed here had produced in 
the plane of temperature and chemical potential a second order transition for low chemical 
potentials and one of first order beyond some critical value of the chemical potential (corre- 
sponding to a so-called tricritical point). This feature induces a related discontinuity in the 
ratio of pion to sigma mass when going beyond the first order transition, that is for higher 
densities. As to the determination of the region in the phase diagram where the sigma is 
no longer kinematically able to decay into two pions, we have calculated it quantitatively 
within the model and found that, as expected, it is indeed confined to a small strip near the 
transition line, and this independently of of the type of transition. We have discussed in 
detail the dependence on the thermal parameters and drawn the expected phase diagram 
behaviours. 
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FIGURES 
FIG. 1. Plot of Mct/Mtt versus T for m = 9.5 MeV . The chemical potential is fixed at 

/i = 70MeV, which is below the value ^ = jit = 76MeV at the tricritical point, for ( = 0.1. 

FIG. 2. Plot of M^/M.,^ versus T for m = 9.5 MeV. The chemical potential is fixed at 
/i = 120MeV , which is above the value fi = jit = 76MeV at the tricritical point, for ( = 0.1. 
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